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The most general solutions of one-dimensional Brownian motion, governed by the 
Chapman-Kolmogoroff functional equation, in the form of series of products of two 
Laguerre polynomials are described. Some analytical and probabilistic properties of these 
solutions are discussed and it is shown that these solutions are different from those given by 
A. N. Kolmogoroff. The limiting behavior of the probability function is investigated. 
A derivation of the corresponding parabolic partial differential equation is obtained. 

1. Introduction 

The mathematical study of Brownian motion and diffusion phenomena in the case of a 
nonuniform fluid led S. Chapman [1] * to introduce the functional equation 

f{p,r+r Q ,r ) t)^{s{p , ,r ) r } t)S{p--p f ) T, i T+p f ) t+r)dp f . (1) 

B. Hostinsky [11] in studying the functional equation of M. Smoluchowski [19] for the 
same prohlem, observed that the functional equation of Chapman was a generalization of the 
functional equation of Smoluchowski, and found the general solution of the Chapman equation 
by using the method of integration by substitutions due to V. Volterra. 

The same author [11] puts the functional equation of Chapman in the more symmetric 
form 

f(XijXz,tijtz)= I f{^l,X2,ti,t2)f(X2yX 3 ,t2,t 3 )dX 2 (2) 

J a 

which is equivalent to that of Chapman for the case where the space in which the Brownian 
motion or diffusion phenomena take place has only one dimension. 

On the other hand, A. N. Kolmogoroff [17] has considered directly a very general func- 
tional equation valid for a domain of any number of dimensions, which reduces to eq (2) in 
the case of a single dimension and which can be considered as a very general form of one of 
the conditions which occurs in the problem of chain probabilities studied by Markoff. 

Moreover, A. Kolmogoroff has given particular solutions of eq (2) and has shown that 
very wide classes of his solutions satisfy linear partial differential equations of the second 
order and parabolic type, and generalized the functional equation of Chapman from the special 
problems of Brownian motion and diffusion phenomena to chain probabilities. The eq (2) is 
commonly called the Chapman-Kolmogoroff functional equation. 

Nevertheless, the most general solution of the Chapman-Kolmogoroff equation has not 
been given. M. Frechet, however, has given very general solutions having forms different 
from the forms given by the preceding authors. 

M. Frechet has found two methods: the first [6] is based on a procedure used in mathe- 
matical physics and the second [7] on reducing the study of the most general solutions which 
are doubly square integrable to the corresponding study in the discontinuous case. The two 
methods provide two distinct forms whose equivalence is neither obvious nor certain. 

The most general solutions of the Chapman-Kolmogoroff equation in one-dimensional 
space have been given by M. Frechet and the present author. In fact, the case of the finite 
interval (0,27r) has been treated by M. Frechet [6], and the case of a region infinite in both 
directions has been investigated in detail [8, 10]. 

i Figures in brackets indicate the literature references at the end of this paper. 
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In this paper some analytical and probabilistic properties of the most general solutions 
of one-dimensional Brownian motion in the interval (0,oo) are described. The corresponding 
parabolic partial differential equation, which constitutes the fundamental equation of Brownian 
motion and diffusion theory, is obtained, and the limiting behavior of the solutions at infinity 
is analyzed. 

2. Statement of the Problem and Formulation of the Fundamental Results 

Let us consider a free spherical Brownian particle (i.e., in the absence of an external field 
force) depending on a finite number v of parameters and which can assume all states Q, forming 
a certain set V (which one can consider as a region of ^-dimensional Euclidean space). One 
assumes then that there exists a probability oo(M,s; v,t) that the Brownian particle passes 
from state M at the instant s to any one of the states of a set v at the later instant t . 

The probability function co(M,s; v,t) is, by virtue of the total probability theorem, an 
additive function of the set v. It is clear that 



u(M,s; v,s)=6(M,v) 

where 

Mev 



fl if, 

) \o if. 



When s<7, one can consider the simple case where the additive function of the set v can 
be represented by a Lebesgue integral 

<*(M,8;r,t)=Jf(M,8;Q,t)dQ (3) 

where dQ is an element of the ^-dimensional space and is a *>-tuple integral. 

When this assumption is made, one finds, in applying total and compound probability 
theorems that the Brownian displacement function /, which govern the ^-dimensional Brownian 
motion, 2 satisfies the Chapman-Kolmogoroff functional equation 

f(M,s; P,t)=f v f(M,s; Q,u)f(Q,u; P,t)dQ (4) 

for s<^u<^t. V denotes a ^-dimensional Euclidean space and /is Lebesgue-measurable in V. 

In the problem of chain probabilities, f(M,s; P,t)dP is the elementary probability that the 
state of a physical system, fixed by v parameters, and represented by M at the instant s } be- 
comes at the later instant t one of the states of a certain elementary set dP of states. However, 
in this interpretation the transition probability or density probability function J(M,s ) P,t) 
must satisfy the probabilistic conditions 



/, 



f(M,s;P,t)>0, (5) 

f(M,s; P,t)dP=L (6) 



If the physical system is moving continuously it has to travel the finite distance MP during 
a finite time t different from s. Then if the states M and P are two different states we have 

lim/(li>; P,0 = 0, M^P. (7) 

If the state P coincides with the state M, the transition probability/ is no more continuous 
and it tends to infinity as t— >s 

lim f(M,s;P,t)= co. (8) 

2 By a ^-dimensional Brownian motion we understand a *>-dimensional strong Markoff process whose Droiections are all one-dimensional 
Brownian motion. 
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Such properties of the transition probability function /, when t— >s, resemble those of the Dirac 
5-function which is a kind of improper function. It follows that the condition (8) must be 
excluded from chain probabilities. Therefore the assumption of the existence of an elemen- 
tary probability f(M,s; P,t)dP cannot be accepted completely and an exception must be 
made for the case s=t, M=P. If /is interpreted as a generalized function, in the sense of the 
distribution theory, then the case s=t and M=P is covered by the statement that / has 
5-function behavior in this case. 

Using the additive theory, both N. Wiener [23] and P. Levy [18] obtained very important 
results concerning the transition probability function. M. Kac [12], S. Kakutani [14, 15] and 
others in a series of papers obtained properties of the Brownian motion in a ^-dimensional 
space and in connection with the Riemann surfaces. An interesting probabilistic discussion 
of the general diffusion process in one dimension is also given by W. Feller [5]. Recently 
M. Kac [13], relying on the properties of Wiener measure in the space of continuous functions, 
showed the disadvantage of Wiener measure theory applied to Brownian motion. 

On the purely physical side the two papers by Uhlenbeck and Ornstein [21] and Uhlenbeck 
and M. C. Wang [22] should be mentioned. 

3. Method of Solution 

To obtain a very general solution of the functional eq (4), M. Frechet [G] tried a method 
of separation of variables in the form 

f(M,s;P,t)=± A t (M,s)B t (P,t) (9) 

and found that the necessary and sufficient condition that the expression (9) should be a solu- 
tion of eq (4) is that, for the values of u such that $</&<£, the functions A t (Q,u) and Bi(Q,u) 
form a biorthonormal system over V, i.e., 

J v A t (Q,u)Bj(Q,u)dQ=S ti where «<r={J Nowise. (10) 

Nevertheless, Frechet realized that the expression (9) does not satisfy the conditions of 
chain probabilities and introduced instead the infinite series 

f(M,s;P,t)=± A n (M,s)B n (P,t) (11) 

n=0 

and proved that the necessary and sufficient condition that the transition probability function 
J(MjS; P,t), satisfying the Chapman-Kolmogoroff eq (4), could be written, in at least one way, 
in the form (11) where the convergence is uniform for s<^t, is that the functions Ai(M,s), 
Bi(P,t) form a complete orthonormal set over V. 

The series solution (11) is supposed to be uniformly convergent for s and t fixed such 
that s<^t, and M,P varying arbitrarily over V. 

4. One-Dimensional Case 

Considering the case of one-dimensional space 0>=1), and with the usual notations x, 
y, and z for M, P, and Q, the Chapman-Kolmogoroff eq (4) can be written in the form 

f(x,s;y,t)=\ f(x,s;z,u)f(z,u;y,t)dz, s<u<Ct. (12) 

Frechet's method led us to take for the solution of eq (12) the infinite series 

f(x,8; y,t)=?0,A n (x,s)B n (y,t) 
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where the functions A n , B n form a biorthonormal set over V. 

The case of the finite interval V(0,2ir) has been treated by M. Frechet [6]. The solution 
of eq (12) on (— oo^oo) in the form of a series of products of two Hermite polynomials has been 
given by the author [8, 10], and its analytical and probabilistic properties are discussed in 
detail. 

The most general solution of eq (12) in the interval (0,oo) is also given by the author 
[8, 9] in the form 

f(x,s; y,t)=±W(x)W(y)d»(s,t)=exp (-^)(xy)l±K 1 L^(xWV(y)e"(s,t),a>-l, (13) 

n=0 \ & / n=0 

where the infinite sequence of Z 2 -f unctions < ^ ( " } (x) > 

fp (x)=\~l exp r-|"pzc-) (x) 

form a complete orthonormal set over (0,eo), and i ( " } (x) is the nth associated Laguerre poly- 
nomial defined by 

r -<x p x fin 

In addition, 

0(M)=^|<i, iovs<t, 

a(s) is a positive increasing continuous function ^0, chosen arbitrarily, and 

_I>+a+l) 
n r(n+l) ' 

where T denotes the gamma function. 

It has been shown [8, 9] that, for s,t fixed such that s<^ and x } y varying arbitrarily over 
(0,°o ), the series solution (13) is absolutely convergent in (0,oo) and uniformly convergent in 
(e,oo ) y e being an arbitrarily small but fixed positive number. 

Using Hardy's formula 

pK l L^{x)L^{y)e-={\~d)-Kexy)-l exp [-^^] /. [^?^} 

where I a denotes the modified Bessel function of the first kind, as the generating function 
for the series solution (13), we get 

/(,,„,«^(l-,,-.exp[-<i|^]/.[ 2 -(M^] (14) 

The modified Bessel function I a (z) is defined by 

oo /l«\a+2» 



po n\(a-\-n)\ 



where z is positive and a is an integer. 

For a a positive integer or zero, a condition which is not incompatible with the assumption 
a> — 1, I a {z) is positive. Thus, according to (14), the transition probability / is always positive 
for any values of x and y in the interval (0,<»). Therefore, the solution (13) satisfies the 
condition (5) provided a takes positive integral or zero values only. 

It has been shown [8, 9] that, as a consequence of Kogbetliantz's inequalities [16], the solution 
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(13) satisfies the special probabilistic condition 

limf(x,s;y,t)=0, 

t->8 

for any finite and positive values of x and y, and also the condition 

lim f(x,s;y,t)= co. 

5. Limiting Behavior of the Probability 

Let us investigate the asymptotic behavior of the probability u(x,s; V,t), defined by (3), 
when the time t tends to infinity. To do so the existence of the limit of the transition proba- 
bility/ defined by the uniformly convergent series 

/(*,«; y,t)= exp [_*±ffj(xy)f { T~\a+l) + ± X^ 1 [j|]] Z*"> (x) Lf > (y) \ (13') 

should be sought. Having in mind that in the factor 0(Sji)=a(s)/a(t) the function a(t) is a 
positive increasing continuous function of t, two cases may occur: 

1. rima(£) = °° 

t->co 

hence 

lim0(s,*)=O. 

Therefore, for x, y, and s fixed and t tending to infinity we have 

limf(x ) s;y ) t) = T- 1 (a+l)(xy)l exp [-(x+y)/2] 

which depends on x and y. It follows from (3) that the limit of the probability (a(z,s; V,t), where 
V is the interval (0,°o), is given by 

lim cofos; V,t) =2*~ +1 r (|+ 1) T" 1 («+ l)al ex P [-f] 

provided a> — 2. This limit depends on the initial state x which is called a nonoscillatory 
case. This terminology is due to M. Frechet, and it means that the probability co tends to a 
limit, in the ordinary sense, when t tends to infinity. 

By assumption, a takes positive integers or zero values, and therefore the condition a> — 2 
is fulfilled. 

2. The only alternative to the case 1 is that 

lim a(t)=d 



where d is finite and positive. In this case 



lim d(s,t)=a(s)ld=ai(s) 

Z->oo 



where 0<ai(s)<l. 
It follows that 



\im J(x,s; y,t)=(xy)l exp f-^1 ± X^&hWW (*)Z?> (y), 
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hence 



Urn u(x,s;V,t)=( (ay)f exp f-^l 2 VkWW (x)L£\y)dy 



where 



n = 



K n =J%g»(y)dy. 



Therefore the limit of the probability co depends completely on the initial state x at the instant s. 
This second case is also a nonoscillatory case. 

The above results can be stated as follows : 
Theorem. If the transition probability f(x,s;y,t) satisfying the Chapman-Kolmogoroff e([ (12), 
is in the form (13), then the limit of the probability co, when t tends to infinity, exists and depends 
on the initial state, i.e., the probability is nonoscillatory, provided a takes zero or positive integral 
values. 

6. Stochastic Partial Differential Equations 

A. Einstein [3] showed that the probability density /satisfies a partial differential equation 
of diffusion type whose coefficient is a certain diffusion constant related to the viscosity and 
temperature of the medium in which Brownian motion takes place, and to the radius of the 
Brownian particle (assumed spherical). 

The existence of the probability density j has been established by many authors from 
different viewpoints, especially by A. N. Kolmogoroff [17] and W. Feller [4], Kolmogoroff 
proved that any function which satisfies the Chapman-Kolmogoroff equation satisfies also a 
second-order partial differential equation of parabolic type, and he obtained a class of particular 
solutions and showed that these, subject to certain special conditions concerning the order of 
the magnitude of moments, satisfy, in variables x and s, the equation 

*<*,.) ^+A(x,) g+g=0 (15) 



where 



B*(x,t)=*1im± f 4 ° (y-x)y(x,t,t+At,y)dy, 

At-^0 *&t J_ ro 

A(x,t)=\im-— (y-x)f(x,t,t+At,y)dy, 

A*->0 A t J _ m 

under the assumption that the ratio of the moment of order 3 to that of order 2 tends to zero 
at the limit At=0. 

Kolmogoroff 's partial differential equation was derived by C. M. Tchen [20] independently 
of Kolmogoroff 's method. 

It can be shown that the solution (13) is different from those discussed by Kolmogoroff. 

To see this, note that if the solution (13) belongs to those of Kolmogoroff, then there would 
exist two functions A and B such that the solution (13) satisfies a partial differential equation 
of type (15). Calculating the derivatives of/, given by (13), we get 

g=(^)f exp [-*+*] g n £ | i- (x)if > (y) 

<^ex P [-^]g£[L;'-(4+(^,)^>(,)+(f+^-^)irw]ir( y ), 



dx 2 ' 



where prime denotes the derivative with respect to x. 
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For x and y varying in (0,<» ) ,t fixed, and s lying in a bounded interval such that s</, the series, 
obtained by differentiating (13) term by term, are uniformly convergent. Substituting the 
above derivatives in (15), one gets 

exp T-^l (xy)f ± ~ [U n (x,s,t)] Z<«> (y) =0 (16) 

where U n (x,s,t) denotes 

Since the relation (16) must hold for any value of y in (0,°°), the expression in the square 
brackets must be zero for any positive integer value of n; that is, after dividing by L^ (x), 



nbd 
6 bs 



^ um<&-t)]^ m§<H §^+(?+^-^)]=»- <"> 



1 7\n 

As - y" is independent of n and x, it is clear that the relation (17) does not hold unless the 

two coefficients A(x,s) and B 2 (x,s) are botli zero. Therefore 6(s,t)=a(s)/a(t) will be inde- 
pendent of s which is contrary to the assumption that a(s) is a positive increasing continuous 
function of s. This proves the statement. 

It can also be shown by computation that the relation (17) is possible if the two functions 
A(x,s) and B(x,s) are both zero. In fact for n=0, we have 

A = x'+a 2 -2a-2x 2 (a+l) 

B 2 2x(x 2 -a) [ b) 

1 

two relations obtained, one gets 



1 y\n 

On the other hand, setting n=l and n=2 in relation (17) and eliminating - r- between the 



A_C 
B 2 ~D 

where 

r= 2(q— x 2 ) x*+a 2 —2a q+1 2x+2(a-x 2 )(x-a— 2) 

x(x-a-iy 4z 2 ' 2 x 3 -2(a+2)x 2 +(a+2)(a+l)x 

and 

2(x— a— 2) 2 . x 



(19) 



D= 



x 2 -2(a+2)x+(a+2) («+l) (x-a-iy2 



The two ratios (18) and (19) are different, hence the relation (17) does not hold unless the 
two functions A and B are both zero. 

Now we propose to obtain a derivation of the corresponding parabolic partial differential 
equation or the so-called backward equation of diffusion theory. For this purpose we can 
use the generating function (14) and get 

I a (z)=dl (1-0) exp [ (1 t(i^ y) ]/(^ ».*) ( 2 °) 

where 

s=2(l -0)- 1 (toy) 1/2 , 
and I a (z) satisfies 

z 2 r:(z)+zU{z)-{a 2 +z 2 )I a (z) = (21) 

prime denotes the derivative with respect to z. 

257 



Taking the derivatives of (20) with respect to x and computing the different terms of (21), 
one gets 

(.• +2 .)7„_ K =<i-,) «p[ (1 tg!! t t y) ] J(1 7i-^ 

Substituting in (21), and dividing by fll (1-0) exp 1" ^^^ y 1 we ha 



lave 



4x2 ^+^ x 



dx 2 
The relations 



f x(l+d) Idf x 2 (l+d) 2 +2x(l-9 2 )-a 2 (l-d) 2 -4dxy 

L 1-0 + J^ f (i-ey ;_u - [ZZ) 



dl a T , dz i dl« T , dz 
-jr~ =±a y- and —-=I a — - 

Oj: ox o# 50 



lead us to express/ in terms of . - and ~- as follows: 

/= 2(1-0) / bf d/\ 

J (l+d)x-a(l-d) 2 +26(l-d)-2d(x+y) \ 50 dxj' 

Substituting the above expression of/ in (22) and setting 

t3(x,s)=2x 2 (l-d)l(l+d)x-a(l-d) 2 +2d(l-6)-2d(x+y))[x 2 (l+d) 2 

+2x(l-d 2 )-a 2 (l-6) 2 -4:dxy]- 1 
and 

y(x,s)=x{2[x(l+e) + (l-d)][(l+d)x-a(l-d) 2 +2d(l-d)-26(x+y)][x 2 (l+e) 2 

+2x(l-d 2 )-a 2 (l-d) 2 -4:dxy]- 1 -l}, 



one gets 

On the other hand 

so that 



0(Mg+7(^)g+*g=O. (23) 

bf__bfb6_a'(s) bf 
bs bd bs a (t) bd 

bf = a(s) a(t) bf = a(s) bf 
bd a(t)a'(s) bs a'(s) bs 



where prime denotes the derivative with respect to s. Hence equation (23) becomes 

^{ +y |/ + ^)|/ = 0. (24) 

dx 2 dx a (s) ds K J 

This derivation presupposes certain differentiability conditions on f and a (s), namely: 
/ must have continuous partial derivatives of order <2 with respect to x and s, and the right- 
handed derivative a'(s) must exist. Therefore the following theorem can be stated: 
Theorem. If the transition probability function j, defined by (13), is of the class C 2 , and the 
right-handed derivative a'(s) of the increasing continuous junction a(s) exists, and is everywhere 
positive and different from zero, then the transition probability f satisfies, in the variables x and 
s, the parabolic equation (24) . 
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It can be observed that the parabolic eq (24) to which the solution (13) satisfies is different 
from equation (15) given by Kolmogoroff. 

The coefficient of ~ in (24) suggests setting a(s)=e 8 , and then pulling a(t)=e\ we have 

OS 

The special case a(s)=e 8 leads us to a particular solution of eq (12) in the form 
f(x,s; y,0=exp [~(x+y)/2](xy)t ± K 1 exp [n(s-t)]L£\x)Lf\y) 

and the eq (24) will be reduced to the normal form 

If we now set, more generally, 

where the variable S increases with s, then 

a(t)=e T 

and 

where /w(<S0 is an increasing continuous function of S. We I hen have 

e=e s ~ T . 

The function /(#,s ; y,t) becomes F(x,S;y,T), and the Chapman-Kolmogoroflf equation will 
still be satisfied if the instants s, t are replaced by numbers S and 7' and I he fund ion F(x,S', y,T) 
can be written as follows: 

F(x,S; y,T)=T F(x,S; z,U)F{z,U; y,T)dz, S<U<T. (27) 

The series solution (13) now becomes 

f(x,s; y,t)=f[x,KS); yXT)]=F(x,S; y,T) 

= exp [-^jp] (ay)* ± K l exp [n(S-T)]LP(x)L™(y). (28) 

When the time t-><x> } T tends either to infinity or (T being negative as well as S) to zero, 
then 6 tends either to zero or e s , so that the limiting behavior of the probability function to, 
defined by 

f* 00 

Jo F ^ ,S] y,T)(hjt 
will still be nonoscillatory. 

There is nothing arbitrary in F(x,S;y,T), defined by (28), which is a particular solution 
of (27). Therefore the most general solution will be in the form 

f(z,8;y,t)=nxM*);V,l(t)] 

where i;(s)=ln a(s) is an arbitrary function varying in (— »,»). 
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